Influence of turbulence on the dynamo threshold 
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We use direct and stochastic numerical simulations of the magnetohydrodynamic equations to ex- 
plore the influence of turbulence on the dynamo threshold. In the spirit of the Kraichnan-Kazantsev 
model, we model the turbulence by a noise, with given amplitude, injection scale and correlation 
time. The addition of a stochastic noise to the mean velocity significantly alters the dynamo thresh- 
old. When the noise is at small (resp. large) scale, the dynamo threshold is decreased (resp. 
increased). For a large scale noise, a finite noise correlation time reinforces this effect. 

PACS numbers: 47.27.Eq, 47.27.Sd, 47.65. +a, 91.25.Cw 
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The process of magnetic field generation through the 
movement of an electrically conducting medium is called 
a dynamo. When this medium is a fluid, the instability 
results from a competition between magnetic field ampli- 
fication via stretching and folding, and damping through 
magnetic diffusion. This is quantified by the magnetic 
Reynolds number Rm, which must exceed some criti- 
cal value Rm c for the instability to operate. Despite 
their obvious relevance in natural objects, such as stars, 
planets or galaxies, dynamos are not so easy to study 
or model. Computer resources limit the numerical study 
of dynamos to a range of either small Reynolds num- 
bers Re (laminar dynamo), modest Rm and Re |lj or 
small Pm — Rm/Re using Large Eddy Simulation 
These difficulties explain the recent development of ex- 
periments involving liquid metals, as a way to study the 
dynamo problem at large Reynolds number. In this case, 
the flow has a non-zero mean component and is fully tur- 
bulent. There is, in general, no exact analytical or numer- 
ical predictions regarding the dynamo threshold. How- 
ever, prediction for the mean flow action can be obtained 
in the so-called "kinematic regime" where the magnetic 
field back reaction onto the flow is neglected (see e.g. 

This approximation is very useful when conduct- 
ing optimization of experiments, so as to get the lowest 
threshold for dynamo action based only on the mean flow 
Rml IF H H IB, It led to very good estimate of the 
measured dynamo threshold in the case of experiments 
in constrained geometries [8j, where the instantaneous 
velocity field is very close to its time-average. 

In contrast, unconstrained experiments are char- 
acterized by large velocity fluctuations, allowing the ex- 
ploration of the influence of turbulence onto the mean- 
flow dynamo threshold. Theoretical predictions regard- 
ing this influence are scarce. Small velocity fluctua- 
tions produce little impact on the dynamo threshold 01 • 
Predictions for arbitrary fluctuation amplitudes can be 
reached by considering the turbulent dynamo as an in- 
stability (driven by the mean flow) in the presence of a 
multiplicative noise (turbulent fluctuations) [ill ]. In this 



context, fluctuations favor or impede the magnetic field 
growth depending on their intensity or correlation time. 
This observation is confirmed by recent numerical simu- 
lations of simple periodic flows with non-zero mean flow 
[lil IfH ] showing that turbulence increases the dynamo 
threshold. 

In the sequel we use direct and stochastic numerical 
simulation of the magnetohydrodynamic (MHD) equa- 
tions to explore a possible explanation, linked with the 
existence of non-stationarity of the largest scales. We 
found that the addition of a stochastic noise to the mean 
velocity could significantly alter the dynamo threshold. 
When the noise is at small scale, the dynamo threshold 
is decreased, while it is increased for a large scale noise. 
In the latter case, the noise correlation time plays a role, 
and reinforces this effect, as soon as it is larger than the 
mean eddy turnover time. When interpreted within the 
Kraichnan-Kazantsev model of MHD flow, these results 
predict that large scale (resp. small scale) turbulence 
inhibits (resp. favors) dynamo action. 

The MHD equations for incompressible fluids are : 



9 t u + u- Vu = -VP + vV 2 u+jxB 
d t B = V x (u x B) + r?V 2 B. 
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Here, u is the velocity, B is the Alfven velocity, P the 
pressure, v the viscosity, 77 the magnetic diffusivity, j = 
Vx (B), and v TG = (sin x cos y cos z, — cosx sin y cos z, 0) 
is the Taylor-Green vortex and f(t) is set by the condi- 
tion that the (1,1,1) Fourier components of the velocity 
remains equal to v TG . The equations are integrated on 
a triply periodic cubic domain using a pseudo-spectral 
method. The aliasing is removed by setting the solution 
of the 1/3 largest modes to zero. The time marching is 
done using a second-order finite difference scheme. An 
Adams-Bashforth scheme is used for the nonlinear terms 
while the dissipative terms are integrated exactly. The 
two control parameters are the Reynolds number Re — 
Vrmshnt/v and the magnetic Reynolds number Rm = 

VrrasUnth, where V rms = (l/3)V2i? = (1/3)V< U 2 > IS 

the (spatial) r.m.s. velocity based on the total kinetic 
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energy E = J E(k)dk and l mt = {Zir/£)E/ J kE(k)dk 
is the integral scale of the turbulent flow. Both v rms 
and li n t fluctuate with time. Thus, viscosity and dif- 
fusivity are dynamically monitored so as to keep Re 
and Rm constant. We have checked that Re is a sim- 
ple linear function of a non-dynamical Reynolds num- 
ber Re cxp — v max Tr/v (usually used in experiments) 
based on maximum velocity and half the simulation box: 
Re = 7.41i?e GX p- In the sequel < X > (resp. X) refers 
to spatial (resp. time) average of X. 

We ran typically four types of simulations : i) DNS- 
MHD, where the full set of equation 0} is integrated 
at 5 < Re < 100 and 5 < Rm < 50 using resolutions 
up to 256 3 ; ii) LES-MHD, where the Lesieur-Chollet 
model of turbulence is used for the velocity equation IB 
a), allowing to explore a case out-of-reach of DNS |12fl . 
namely Re = 500, 5 < Rm < 100; iii) kinematic sim- 
ulations ; iv) kinematic-stochastic (KS) simulations. In 
these last two cases, only the induction equation Jfl-b) 
is integrated with u set to a given velocity field. In the 
kinematic case, it corresponds to the stationary velocity 
field u(i?e) obtained through time-average of a stable so- 
lution of the Navier-Stokes equations with Taylor-Green 
forcing, at fixed Reynolds number. This procedure is 
complicated by the presence of hydrodynamic instabili- 
ties at low Reynolds number |l4| , which impose very long 
simulation time (typically over 1000 s, i.e. 400 eddy turn 
over times at Re = 46) to ensure convergence towards an 
asymptotically stable solution. The average is then per- 
formed over several (typically 200) eddy-turnover times. 
In the KS case, the velocity field u = u(Re) + v'(ki,T c ) 
is the sum of a time averaged velocity field at a given Re 
and of an external Markovian Gaussian noise, with fixed 
amplitude v' , correlation time r c and typical scale ki. In 
both kinematic simulations, the magnetic Reynolds num- 
ber Rm is computed by using the rms velocity and inte- 
gral scale of u. In the deterministic case, this amounts 
to use V rm s = (l/3)\/ < u 2 > and Li nt the (spatial) rms 
velocity and integral scale of the time-averaged velocity 
field, therefore respecting the experimental procedure fol- 
lowed in optimization of dynamo experiments 0,0, 0,0- 

For each type of simulation, we fix Re ( v', t c and ki, if 
needed), vary Rm and monitor the time behavior of the 
magnetic energy < B 2 > and the finite-time Lyapunov 
exponent A = 0.5 d t < ln(_B 2 ) >, where the average is 
taken over the spatial domain. Three types of behaviors 
are typically observed fl4| : i) no dynamo : the magnetic 
energy decays, the Lyapunov converges towards a finite 
negative value; ii) undecided state : the magnetic en- 
ergy remains at a low level, with intermittent bursts of 
magnetic energy and oscillation of the Lyapunov, so 
that no fit of the Lyapunov exponent can be obtained; 
iii) turbulent dynamo : the magnetic energy grows with 
positive Lyapunov, and, in the DNS-MHD or LES-MHD, 
reaches a nonlinear saturated regime. 

From the values of the Lyapunov in the turbulent dy- 



namo and no dynamo regime, one may derive the crit- 
ical magnetic Reynolds number Rm c (Re), solution of 
A(i?e, Rm c ) = 0, through a standard interpolation pro- 
cedure. 
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FIG. 1: Simulation parameter space. Square refer to DNS- 
MHD and LES-MHD simulations, and shaded areas to win- 
dows of dynamo action for kinematic simulations with mean 
flow. □ : no-dynamo case ; EB : undecided state; ■ : dynamo 
case ; — Rrn t c urb ; Rm^ IF ; — • — ■ end of the first kine- 
matic dynamo window ; • ■ ■ beginning of the second kinematic 
dynamo window. Shaded areas indicate explored windows of 
dynamo action for the mean flow. 



A summary of our exploration of the parameter space 
is provided in Fig. 1, for the non-stochastic simulations, 
where the only control parameters are Rm and Re. We 
did not detect any dynamo at Re = 2. Between Re = 4 
and Re = 6, we observed heterocline dynamos, oscillating 
between a non-dynamo and a dynamo state. The win- 
dow 2 < Re < 4 has been studied in 0, where decreas- 
ing critical magnetic Reynolds number has been found. 
For 4 < Re < 100, we found that the critical magnetic 
Reynolds number for dynamo action in a turbulent flow 
i?m'" rf> increases with the Reynolds number, in quanti- 
tative agreement with the result obtained in the same 
geometry, but with a different forcing (at constant force 
instead of constant velocity) 0] . Our LES-MHD simu- 
lation confirms the saturation of the dynamo threshold 
at large Reynolds number already observed in constant 
force simulations 0. For the mean flow, we have actu- 
ally detected at least two windows of dynamo actions : 
one, independent of Re, starting above Rm^ IF « 6 and 
centered around Rm =10, with real Lyapunov (station- 
ary dynamo); a second, occurring at larger Rm, varying 
with Re, with complex Lyapunov (oscillatory dynamo). 
One sees that i?m*" rb varies across these two windows 
and always exceeds Rm^ 1F . In the sequel, we show that 
the increase and saturation of i?TO* Mrb is not due to a 
crossing between the two dynamo modes, but to the in- 
fluence of non-stationary large scales over Rm^f F . 

To make an easier connection between DNS and KS 
simulations, we introduce a parameter that quantifies the 
noise intensity, S — < u 2 >/ <u 2 >. This parameter de- 
pends on the noise amplitude, as well as its correlation 
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FIG. 2: a) Noise intensity 8 = < u 2 >/ < u >, as a function 
of the Reynolds number, measured in our DNS simulations 
under the dynamo threshold, b) Energy spectrum of the ve- 
locity field in the DNS, at Re — 46, for different average 
period T: dotted line : T = ; short dashed line : T = 75 s ; 
long-dashed line : T = 150 s ; continuous line : T = 300 s. 




time and characteristic scale, and need to be computed 
for each stochastic simulation. It can also be computed 
in the direct simulations, and is found to depend onto 
the Reynolds number, increasing from a value of 1 at low 
Reynolds number, until about 3 at the largest available 
Reynolds number (Figure |2J-a) . Note that 8 — 1 is just 
the ratio of the kinetic energy of fluctuations onto the 
kinetic energy of the mean flow. In the sequel, the com- 
parison between the KS and DNS-MHD simulations will 
therefore be made using 8 as the control parameter. An- 
other interesting information can be obtained from the 
energy spectrum of the velocity field, as one averages over 
longer and longer time scales (Figure |2J-b). One sees that 
during the first period of average (typically, a few eddy 
turn-over time, i.e. about 5 to 10 s), one mainly removes 
the fluctuations at largest scales, while the remaining av- 
erage mostly removes small scales (over time scales of the 
order of 50 to 100 eddy-turnover times, i.e. 300 s). 

In the sequel, we explore the influence of both type 
of fluctuations through the KS simulations, by consider- 
ing noise at large (fcj = 1) and small scale (k] — 16), 
with correlation time ranging from to 50 s. Since the 
kinematic dynamo threshold is essentially constant for 
all values of the Reynolds number we explored, we first 
focus on the study of the case where the time-averaged 
field is fixed as u(Re = 6) and vary the noise ampli- 
tude, characteristic scale or correlation time, to explore 
their influence on the dynamo threshold. An example 
of our exploration of the parameter space is provided 
in Figure |3| for different kinds of noise and u(Re = 6). 
Note that by using our external noise, we are able to 
produce noise intensities of the order of noise intensities 
measured in experiments (8 ~ 10 at Re ~ 10 6 for the 
von Karman flow), and that are out of reach of DNS. 
For low correlation time or injection scale, we are actu- 
ally able to follow the deformation of the two windows of 
dynamo action. One sees that a noise does not destroy 
them, but rather distorts them. In the case where the 
noise is at small scale {kj = 16), the two windows are 



FIG. 3: Parameter space for noise at Re — 6 for different 
noise parameters : a) r c = 0, ki = 1; b) r c = 0.1 sec, ki — 1; 
c) r c — 0, ki — 16; d) r c = 0.1 sec, ki — 16. □ : no-dynamo 
case; EB : undecided state; ■ : dynamo case. The full lines 
are zero-Lyapunov lines. 



slightly tilted downwards, while they are lifted upwards 
in the case of large scale noise (kj = 1). The influence of 
the noise onto the first dynamo bifurcation (the dynamo 
threshold) can be summarized by plotting the critical 
magnetic Reynolds numbers as a function of the noise 
intensity (Fig. EJ-a). Large scale (resp. small-scale) noise 
tends to increase (resp. decrease) the dynamo thresh- 
old. Furthermore, one sees that for small scale noise, the 
decrease in the dynamo threshold is almost independent 
of the noise correlation time r c , while for the large scale 
noise, the increase is proportional to r c at small r c . At 
t c > 1 sec-one third of the mean eddy-turnover time-, all 
curves Rm c {8) collapse onto the same curve. We have 
further investigated this behavior to understand its ori- 
gin. Increasing 8 first increases of the flow turbulent 
viscosity v rms li nt with respect to its mean flow value 
V rms Lint- This effect can be corrected by considering 
Rm* — Rm c V rms Li nt / v rm Ji nt . Second, an increase of 
8 produces an increase of the fluctuations of kinetic en- 
ergy, quantified by S 2 = \J < u 2 > 2 — < u 2 > 2 /< u 2 >. 
This last effect is more pronounced at kj = 1 than at 
fcj = 16. It is amplified through increasing noise correla- 
tion time. We thus re-analyzed our data by plotting Rm* 
as a function of 82 (Fig. IH-b). All results tend to collapse 
onto a single curve, independently of the noise injection 
scale and correlation time. This curve tends to a con- 
stant equal to Rm^ IF at low 82- This means that the 
magnetic diffusivity needed to achieved dynamo action 
in the mean flow is not affected by spatial velocity fluc- 
tuations. This is achieved for small scale noise, or large 
scale noise with small correlation time scale. In contrast, 
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the curve diverges for 82 of the order of 0.2, meaning 
that time-fluctuations of the kinetic energy superseding 
20 percent of the total energy annihilate the dynamo. 

We now turn to detailed comparison of dynamo thresh- 
olds obtained in KS simulation with DNS-MHD case. In 
Fig. El we show that the dynamo threshold obtained at 
fcj = 1 , for t c > 1 as a function of the noise intensity cor- 
responds to the DNS-MHD dynamo threshold. Note that 
the noise intensity 8 saturates past a Reynolds number of 
about 100 (Fig. EJa), thereby inducing the saturation of 
the critical magnetic Reynolds number at large Reynolds 
number. To check that our results are not affected by the 
choice of u, we ran additional KS simulations with u com- 
puted at Re = 25,46 and 100. Since the computational 
cost in these cases is much larger than in the case Re — 6, 
we focused on the case where the noise has a correlation 
time r c = 1 or 8 sec. and injection scale fc/ = 1 and 
only computed the critical magnetic Reynolds number for 
the level of noise reached by the DNS at that Reynolds 
number. Figure [S] shows that the dynamo threshold co- 
incides with the dynamo thresholds both of KS simula- 
tions at Re — 6 and of the DNS, indicating that a large 
scale noise is probably responsible from the increase of 
Jl m turb Reynolds number. A physical identification 
of the dynamics of the velocity fluctuations playing the 
role of this noise can be performed by visual inspection of 
the turbulent velocity field. One observes that the large 
scale vortices generated by the Taylor-Green forcing are 
not exactly stationary, but wander slightly with time. A 
similar large-scale non-stationarity has been observed in 
the shear layer of Von Karman flows ^| 0|. This pro- 
cess is approximately reproduced by a large scale noise 
with sufficiently long correlation time, and may therefore 
been held responsible for the significant increase of the 
dynamo threshold. 

Our work suggests that it might not be so easy to 
achieve turbulent dynamos in unconstrained geometries, 
with large scale non-stationarity. In the experiments, a 
necessary ingredient for dynamo action could therefore 
be a monitoring of the large scale, so as to keep them as 
stationary as possible. In geo- and astrophysical flows, 
this role could be played by the Coriolis force. Our work 
also indicates that a well chosen noise can be used in place 
of the actual turbulent velocity fluctuations to compute 
the dynamo threshold, at a much lower computational 
cost. In some sense, a kinematic-stochastic simulation 
can therefore be seen as a turbulent model and might be 
useful in the astro- or geophysical context. 
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FIG. 4: Evolution of the dynamo threshold for KS simulations 
with u(Re — 6). a) Rm c as a function of 5 and b) Rm* as 
a function of 82 for different noise parameters : k = 1 : □ 
t c = ; □ t c — 0.1 sec; B r c = 1 sec; EH t c = 8 sec ; ■ r c = 50 
sec ; k = 16 : o r c = ; t c = 0.1 sec; • r c = 50 sec. 
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FIG. 5: Evolution of the dynamo threshold as a function of 
8 for the DNS (squares, same meaning as in Fig. 1 and 3) 
and KS simulations with ki = 1 using u(Re = 6), r e = 1 
sec (• connected with line), and u at the Reynolds number 
corresponding to the DNS with equivalent 8 (o for r c = 1 sec 
and for r c = 8 sec). 
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